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Abstract 

We introduce the harmonic oscillator on the Lobachevsky plane with the aid of 
the potential V(g) = (a 2 w 2 /4) sinh(^/a) 2 where a is the curvature radius and g 
is the geodesic distance from a fixed center. Thus the potential is rotationally 
symmetric and unbounded likewise as in the Euclidean case. The eigenvalue 
equation leads to the differential equation of spheroidal functions. We provide a 
basic numerical analysis of eigenvalues and eigenfunctions in the case when the 
value of the angular momentum, m, equals 0. 



^ ' 1 Introduction 

On 

This short contribution is dedicated to the memory of Vladimir Geyler, our dear 
colleague and friend who, to our great sorrow, passed recently and very unexpectedly 
away. Among his numerous scientific activities in the framework of mathematical 



C — ! physics one can point out his interest in the influence of nontrivial geometries on the 

properties of quantum mechanical systems which was very intensive especially in the 
last period. Let us mention just a few of his results focused on this sort of problems: the 
^ effect of surface curvature on magnetic moment and persistent currents [6j, the general 

form of on-diagonal singularities of the Green function for Schrodinger operators on 
curved spaces [5], and zero modes in a system of Aharonov-Bohm solenoids on the 
Lobachevsky plane [5]. 

In one of the very last scientific contributions due to Vladimir Geyler 1 1 j . on which 
we had the pleasure to collaborate, a quantum dot with impurity in the Lobachevsky 
plane is discussed. This paper was intended as a step towards an extension to the 
hyperbolic plane of results derived in [I] where a similar problem is treated in the case 
of the Euclidean geometry. The current article aims to complete the discussion from 
[TU] by a short numerical analysis of spectral properties of the harmonic oscillator on 
the Lobachevsky plane. In particular, we are interested in the effect of the surface 
curvature on these properties. 
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It has been observed already some time ago that various solvable quantum me- 
chanical systems admit an explicit treatment also in the case of hyperbolic geometry, 
though with interesting nontrivial modifications (see, for example, [El [I])- Also the 
spectral problem for the quantum harmonic oscillator on the Lobachevsky plane, as 
we introduce it below, leads to a differential equation which is well known from the 
theory of special functions, namely to the differential equation of spheroidal functions. 
It should be emphasized, however, that the history of spheroidal functions is much 
more recent than that of more traditional special functions like Bessel functions or 
Legendre polynomials. For example, one of the basic monographs devoted to these 
functions appeared only in the fifties of the last century [TTJ. Also the notation is 
still not fully uniform. One can compare, for example, [TT] with [2]. Here we follow 
the latter source. Furthermore, there are cases of values of parameters on which the 
spheroidal functions depend that have not been fully investigated. In this connection 
we note that particularly the values of parameters which are of interest for our model 
are treated in the textbooks in a rather marginal way. Also effective numerical algo- 
rithms to evaluate the spheroidal functions seem to be rather tedious to create and 
not available for all cases. 

These circumstances make the numerical and qualitative analysis of the spectrum 
more complicated than one might expect at first glance. The numerical results pre- 
sented in this contribution were derived with the aid of the computer algebra system 
Mathematica. On the contrary to other "more traditional" special functions, the 
spheroidal functions were implemented in this system only in the version which was 
the very last one at the time of this writing. 



2 The Hamiltonian 

Denote by a, < a < oo, the so called curvature radius which is related to the scalar 
curvature by the formula R = — 2/a 2 , and by (g, 0), < g < oo, < < 2tt, the 
geodesic polar coordinates on the Lobachevsky plane. Taking the mass equal 1/2 and 
h = 1, the Hamiltonian has the form 

( d 2 I fo\ d 1 d 2 1 \ 

ydg 2 a \aJ dg a 2 sinh(f) <90 2 4a 2 J 

and is defined in L 2 ((0, oo) x S 1 , a smh(g/ a)dg dcj)) on its natural domain. Here V(g) 
is the potential of the harmonic oscillator on the hyperbolic plane. Its choice is am- 
biguous. Of course, one requires that in the flat limit, a — > oo, V(p) becomes the usual 
potential of the isotropic harmonic oscillator on the Euclidean plane, V°°(g) = u 2 g 2 /4:. 
At least three choices of V(g) have been discussed in the literature so far. 
A distinguished possibility is 

Vi(o) — _ a 2 uo 2 tanhf — 

since in that case one gets a Hamiltonian with a comparatively rich symmetry, a so 
called super-integrable model, this is demonstrated in detail on the level of classical 
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mechanics in |12j . see also |13j . Notice that this potential is bounded and thus the 
spectrum of the corresponding quantum Hamiltonian contains an absolutely contin- 
uous part, the interval [£'o(a),oo[. The spectral analysis for this case is covered by 
paper [6J. It turns out that E (a) ~ a 2 as a — > oo, and so the continuous part disap- 
pears in the flat limit. Moreover, as expected, the number of eigenvalues in the curved 
case is finite but tends to infinity in the flat limit, and the eigenvalues approach the 
discrete spectrum of the flat Hamiltonian. 

Further, a nonlinear quantum model with a potential of the type 

Q 2 

V 2 (g) = const 



l + Xg 2 

depending on a parameter A is discussed in [7]. The treatment of this model also relies 
heavily on symmetry properties which allow for explicit solutions. 

In the current paper we stick, however, to a choice made in and otherwise not 
discussed in the above cited papers. We set 

V(g) = ^a 2 u 2 smh^J . 

With this choice, the potential is unbounded, likewise as in the Euclidean case. Conse- 
quently, the corresponding Hamiltonian H has a purely discrete spectrum depending 
on the parameter a. 

Next we rescale the Hamiltonian letting 

H = a 2 H. 

Substituting £ = cosh(^/a) we obtain 

*-(i-rt$-*» + (i-rt-£ + -£«.-i,-i 

in L 2 ((l, oo) x S 1 , d£d</>). Using the rotational symmetry which amounts to a Fourier 
transform in the variable (f>, H may be decomposed into a direct sum as follows 



H = 2__, H mi 



m=—oo 



For all m G Z, H m acts in L ((1, oo), d£). 



3 Eigenvalues and eigenfunctions 

Denote by H m , m G Z, the restrictions of H to the eigenspace of the angular mo- 
mentum with value m. We prefer to regard H m as an operator in L 2 (W + ,dg). This 
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assumes an obvious unitary transform to pass from the measure asinh(g/a)dg in the 
L 2 space on M.+ to the measure dg. On the general grounds, the spectrum of H m is 
semibounded below, discrete and simple 0, EH] ■ Set 



a 2 uj 



Q 
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If E™, n = 0,l,2,..., is the nthe eigenvalue of H m then 

1 ~ E m 

is the nth eigenvalue of H m . 

The eigenvalue equation H m tp = Eip written in the coordinate £ takes the form 

where 

Differential equation (TjQ) is known in the theory of special functions as the differential 
equation of spheroidal functions [2,[TT]. The notation used below follows the source [2]. 
All parameters in ([T|) are in general complex numbers. There are two solutions that 
behave like £^ times a single-valued function and times a single- valued function 

at oo. The exponent v is a function of A, 9, m, and is called the characteristic exponent. 
Usually, it is more convenient to regard A as a function of u, m and 9. We shall write 
A = A™(#). The functions \™(9) obey the symmetry relations 

W) = K m (8) = XHu-iV) = A^LiW- (2) 

A first group of solutions (the so called radial spheroidal functions) is obtained as 
expansions in series of Bessel functions. They are denoted S™ (£,#), j = 1,2,3,4. 
Several relations between the radial spheroidal functions are known, notably 

S™®(£,9) = Sf 1 \£,9)+iSf 2 \£,9) = -ie-^S^ 3 \^,9) 

and 

S™ {3) (£,9) = . 1 (s m i 1 } l {^9) + ie-^S^ l \^9) 
«cos(i/7r) V / 

Also the asymptotic behavior of radial spheroidal functions is well known. It turns 
out that among these functions the only one which is square integrable near infinity 
is S„ m (£,9). For q > it holds 

(e, 4) = ex K-^ - 4 t ) (1 + 0(rl)) 

as £ — ► +oo. 
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To control the behavior near the singular point £ = 1 one uses expansions in series 
in Legendre functions. These solutions are called angular spheroidal functions and are 
denoted Ps™(£,9) and Qs™(£,,9). There exist relations between these two types of 
spheroidal functions. In particular, it holds 

where V™(6) is a coefficient independent of £. Consequently, 

sr (3) (£,0) = . 7 s {vZ-MQs^^,e)+ie- i ^v v m (e)QsZ- 1 ^,e)) ■ 

ZCOS(z/7T) 

The angular spheroidal functions asymptotically behave, as £ — > +00, in the fol- 
lowing way: 

"> = r( i-m)cW ( (1 + ° K " 1}) 

where s™(6) is some coefficient (for its definition we refer to [2]). Particularly, we have 

1 



Furthermore, 



and 



It follows that 



Q«;K,9) = - s Ljiog«-i) + o(i) 



"» = " W (%T + 1oeK " j) + 0(1) - 

We conclude that = ,C (3) (£, -g 2 /4) is an eigenfunction and, at the same 

time, E™ = — A™(— g 2 /4) — 1/4 is an eigenvalue of H m if and only if t/>(£) has no 
singularity as £ — > 1+. Further we restrict ourselves to the case m = 0. The function 
A°(0) is real for v — —1/2 + it, t G R, and becomes negative if t is sufficiently large. 
The eigenvalues of H are the numbers — A u 2+ it{~ <? 2 /4) — 1/4 where t solves the 
equation 

-W2 y -W-g 2 / 4 ) _ ^ jWjg = n 

«°-i/ 2+it (-? 2 /4) - - 1/2 ^(-^ 2 /4) • 

Because of symmetry (J2J) it suffices to consider t > 0. Equivalently, in the numerical 
computations the following form of this equation proved itself to be rather effective: 

gb.G-v ae = - 
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Finally, we present a few numerical results. Figure[T]depicts several first eigenvalues 
of H . More precisely, the figure contains plots of eigenvalues E®(q)/uj = E®(q)/(2q) 
for n — 0,1, 2. These functions of the variable q = a 2 u/2 decrease monotonically, and 
in the limit q —>■ oo they approach the eigenvalues in the flat case, namely 2n + 1. On 
the other hand, for q 0+ the eigenvalues tend to infinity. 

Next we plot the corresponding eigenfunctions for several values of the parameter 
q. We have chosen q = 0.5, 5, oo. The last value corresponds to the flat case. To 
reach this goal it was necessary to fix the value of uj. We have taken u = 1. Then the 
chosen values of q correspond to the curvature radius a 2 = 1, 10, oo. An eigenfunction 
of H equals (£, — g 2 /4) times a normalization factor for a suitable value of 
v G — l/2 + iR+. The plotted eigenfunctions are expressed in the variable g. Thus the 
transformed eigenfunction of Hq in L 2 (M + ,d^) equals 

^ ) = U sinh (!)J n cosh (f))- 

The fact that we stick to the same Hilbert space in all cases, namely to L 2 (R + ,dg), 
facilitates the comparison of eigenfunctions for various values of the curvature radius. 
The same Hilbert space is used also in the limiting Euclidean case when a = oo. The 
eigenfunctions for the flat case are, of course, well known from textbooks. Let us recall 
that in the Euclidean plane we have 

€(^) = (-irv^^(y)exp(-^). 

Looking at figures [21 El and H] depicting the eigenfunctions indexed n — 0, 1 and 2, 
respectively, one observes clearly that a more curved space (corresponding to smaller 
a) forces the eigenfunctions to concentrate more closely to the center of the harmonic 
oscillator. 
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Figure captions 

FIGUREHJ Plots of E° n (q)/uj = E Q n {q)/{2q), n = 0, 1, 2, where E° n is the nth eigenvalue 
of Hq and q = a 2 u/2. 



FIGURE [2j The normalized eigenfunction of H , iJ)q(q) G L 2 (IR + ,d^), for the values 
to — 1 and a 2 = 1 (solid), a 2 = 10 (dashed), a 2 = oo (= the flat case, dotted). 



FIGURE 3. The normalized eigenfunction of H , ipx(g) G L 2 (IR + , df?), for the values 
to = 1 and a 2 = 1 (solid), a 2 = 10 (dashed), a 2 = oo (= the flat case, dotted). 



FIGURE 4. The normalized eigenfunction of H , ip^io) ^ -^ 2 (^+)d£>), for the values 
lo = 1 and a 2 = 1 (solid), a 2 = 10 (dashed) and a 2 = oo (= the flat case, dotted). 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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Figure 4: 
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